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Abstract 

We construct a number of new (v;r, s;X) supplementary difference sets (SDS) with v 
odd and A = (r + s) — (v — l)/2. In particular, these give rise to D-optimal matrices of 
' the four new orders 206, 242, 262, 482, constructed here for the first time. 

(N ' 

1 Introduction 

! Let v be an odd positive integer and consider (—1, l)-matrices H of order 2v. Ehlich's bound 
states that 

Q ' det{H) <2 v (2v-l){v-l) v - 1 . 

D-optimal matrices are 2v x 2v (—1, l)-matrices that attain Ehlich's bound, i.e. they have 
maximal determinant. Ehlich also proved that if A and B are circulant (—1, l)-matrices of 
order v such that 

AA T + BB T = 2(v - l)I v + 2J V (1) 



(where is the v xv identity matrix and J v is the v xv matrix with all elements equal to +1), 



. then the matrix 

^ H=( i T f T ) (2) 



-B 1 A 1 

has maximal determinant. Such A and B can be constructed by using cyclic SDSs, say (X, Y), 
with parameters (v;r,s;X). Thus |X| = r, \Y\ = s and A = r + s — We say that D- 
optimal matrices ()2]) are of circulant type and that such SDSs are D-optimal. By pre- and 
post-multiplying equation (CQ) with J v , one obtains that a 2 + b 2 = 4v — 2, where a and b are row 
sums of A and B, respectively. Each integer solution (a, b) of this Diophantine equation can be 
^ \ used to obtain feasible parameters r and s for the required SDS. Namely, we need a = v — 2r 
and b = v — 2s. By a normalization, we may assume that < a < b which implies that r > s. 
A comprehensive table of all odd v < 100 for which D-optimal SDSs are known, can be found 
in [3]. Additional values of v < 100 for which D-optimal matrices of order 2v have been found 
subsequently, can be found in [12J. There are two infinite series of D-optimal matrices, one for 
v = q 2 + q + 1 where q is a prime power, see [H], and one for v = 2q 2 + 2q + 1 where q is an 
odd prime power, see \17\ . Apart from these two infinite series, the only odd values of v > 100 
for which D-optimal matrices of order 2v are currently known are v = 113, see [3l v — 145, 
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see [3j |8], v = 157, see [8j [TO] and v = 181, see [HJ, [15] . We point out that the cases v = 157 and 
v = 181 were overlooked in |12j . 

In this paper we construct several new D-optimal SDSs for seven new parameter sets: 

(63; 29, 24; 22) (93; 45, 37; 36) (103; 48, 42; 39) 

(103; 46, 43; 38) (121; 55, 51; 46) (131; 61, 55; 51) 
(241; 120, 105; 105). 

In particular, this means that we have constructed the first examples of D-optimal matrices of 
orders 206, 242, 262 and 482. D-optimal matrices of order 126 and 186 constructed previously 
in [2] and [3] were also derived from D-optimal SDS but with different parameters, namely 
(63; 27, 25; 21) and (93; 42, 38; 34), respectively. However the smallest order n = 2 (mod 4) 
for which a D-optimal matrix is still unknown remains n = 138. For a comprehensive list of 
unknown D-optimal SDSs see Section 5. 

We also prove two theoretical results for D-optimal matrices. The first is a generalization of 
the horizontal and vertical constraints for D-optimal matrices, recently proved in [13]. The 
second is the property that if a D-optimal SDS solution is formed by taking union of cosets of a 
certain subgroup of the group of units Z*, then the power spectral density of the corresponding 
sequences must be constant on the cosets. (See below the definition of the power spectral 
density.) This property can be restated as saying that several power spectral density values 
corresponding to the first rows of the circulant matrices A, B in (J2J), must be equal. Both 
our theoretical results can potentially be used as subroutines to discard candidate sequences or 
candidates SDSs, within algorithmic schemes that search for D-optimal matrices. 
Let us recall the definition of the power spectral density (PSD) of a complex sequence X = 
x , xi, . . . , a^-i- Let to = e~^~ be a primitive u-th root of unity and define the discrete Fourier 
transform (DFT) of X by 

DFT x {k) = J2 x ^ ak , k = 0,...,v-l. 

Then the power spectral density values PSDx{k) are defined by 

PSD x (k) = \DFT x {k)\ 2 , fc = 0,l,...,v-l. 
From now on we work only with {±l}-sequences. The well-known symmetry properties 
DFT x (v - k) = DFTx(k) and PSD x (v - k) = PSD x (k), k = 0, 1, . . . , v - 1 
are used within most algorithmic schemes that search for D-optimal matrices. 

2 Generalized horizontal and vertical constraint 

We state and prove a generalized constraint, that yields the horizontal and vertical constraints 
proved in [13] for v = (mod 3), as corollaries. In particular we show that actually any divisor 
d of v gives rise to a Diophantine constraint. 
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Theorem 1 Let the first rows of the circulant matrices A, B in be denoted as [do, • • • , a v ~i] 
and [bo, ■ ■ ■ , &u-i] respectively, where a^, 6j G { — 1, +1}, « = 0, . . . , v — 1. Lei v = dm and set 

A? = a i + a i+ d + ' ' ' + Oj+(m-l)d 

-^i = + fy+d + 1" &j+(m-i)d 

/or j = 0, . . . , d — 1. If A,B define a circulant D-optimal matrix, then 

d-l 



and 



J2(A 2 J +Bf) = 2(v + m-l) (3) 
i=o 

£(A fc A, + B fc B,) = v - m. (4) 



Proof 

By abuse of notation we denote by A, B the circulant matrices and the sequences that are 
specified by their first rows. An argument similar to the one given in [7], shows that the 
hypothesis that A,B define a circulant D-optimal matrix implies that the sum PSDa{s) + 
PSDb(s), s 7^ 0, is equal to the constant 2v — 2, i.e. 

E' 2nks\ f^-k . 2nks\ /\-^, 27rA;s\ / ^-4 . . 2nks\ 
a k cos I + y. a k sm I + bk cos J + y. bk sin J = 2v — 2 

k=0 V J \k=0 V J \k=0 V J \k=0 V J 

for all s G {1, 2, . . . , v — 1}. (For s = the sum is equal to 4t> — 2.) Considering the particular 
values s = mr for r — 1, 2, . . . , d — 1 we obtain: 

27rjbV 27r£rV /^i 27rA:r\ 2 , 27r£r\ 2 

^ ^ a fc cos — — I + I y j a k sm — — I + I ^ j b k cos — — I + I y b k sm — — I = 2u-2. 

vfc=0 / \/c=0 / \fe=0 / \fc=0 / 

(5) 

Noting that 

E 1 27rfcr\ / ^— I 27rA;r\ / ^-i, 2tx jr\ ( i-X 2ixjr 

a fe cos— I + l^ flfc 8in— I = + £ A, sm — 

vfe=o / \fc=o / \i=o / \i=o 



j'=0 k<l 

(and similarly for the B sequence), we see that (J5J) can be rewritten as 

+ Bj) + 2 £(,4*4 + S*S0 cos 27r(/ ~ fc)r = 2« - 2. (6) 
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The hypothesis that A, B define a circulant D-optimal matrix implies that the sum of the 
squares of the sums of the elements of the sequences A and B is equal to the constant 4v — 2, 
i.e. 

'd-i \ 2 /d-i \ 2 

,fc=0 / \k=0 J 

which can be rewritten as 

d-i 

Y,{A 2 j + S?) + 2^2(A k A t + B k Bt) =Av-2. (7) 

j=Q k<l 

By taking the sum of (jSJ) over r = 1, . . . , d — 1 and of ((7j), we obtain the equality (j3J). Now (j4]) 
follows from (J7J). 

□ 

Taking the difference of flT]) and we obtain 



0<fc<Z<<2 

which can be rewritten as 

(A fc A + fifcfiO sin 2 1 - ; = -. (8) 

0<k<l<d 

We remark that assuming v = (mod 3) and applying theorem [1] for c? = 3 and d — | yields 
immediately the vertical and horizontal constraints proved in [T3] . 



Corollary 1 

(a) Assuming v = (mod 3) and taking d = 3 in theorem^ we obtain the vertical constraint: 

A 2 Q + Al + A 2 2 + Bq + B 2 + B 2 = 8m-2 

m—1 m—1 



w/iere = a M+j , Bj = ^ hi+j,3 = 0, 1, 2. 



i=0 i=0 

(b) Assuming v = (mod 3) and taking d = ~ m theorem^ we obtain the horizontal constraint: 



i-i 



J2(A* + B?) = 2v + 4 



i=0 



where = a* + a i+(i + a i+2( i, 5, = b { + b i+d + 6 i+2 d, i = 0, . . . , d - 1. 

Let us also illustrate the interesting phenomenon that in case the irrationalities (that the 
presence of the sine function in (jSJ) potentially entails) do not cancel out, then one obtains 
more than one Diophantine equation. 
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Corollary 2 Assuming v = (mod 5) and taking d = 5 and r = 1 in we obtain the 
constraints: 

( A kA + B k Bi) = 4m and ^ (A k A t + B k B t ) - ^ ( AkAl + = 

0<fc<Z<5 l-k e{2,3} Z-fce{l,4} 

m—1 m—1 

w/iere A,- = ^ a ^+v B i = Z^l b ^+i^ = 0, 1, . . . , 4. 

1=0 «=o 

Proof 

For d — 5, v = 5m, the summation in (jSJ) for r = 1 extends over the 10 possible (A;, /) pairs with 

< k < / < 4 and m each case we have that sin is either equal to sin — = 

~ 5 H 5 8 

2-7T 5 -|- 

or equal to sin 2 — = — - — . Since 1 and v^5 are linearly independent over the rationals, we 
obtain the aforementioned constraints. 

□ 



3 Power spectral density constancy over orbits 

First we introduce some notations and terminology. Let Z v be the ring of integers mod v, i.e 
Z v = {0, 1, . . . ,v — 1}. Let Z* be the group of invertible elements of Z v , i.e. Z* = {k 6 Z v : 
gcd(A;,t>) = 1}. The order of Z* is equal to <f>(v), where denotes the Euler totient function. 
Let H ^ Z* be a subgroup of Z*. Then H acts on Z v and we denote the orbits of this action 
by 

O 1 = {0},O 2 ,...O m . 
Thus we have the disjoint union relationship Z v = 0\ U O2 U • • • U O m . 

The method used in this paper as well as in [21 [31 El El QUI US] constructs solutions for circulant 
type D-optimal matrices by expressing the corresponding SDSs as unions of certain orbits 
associated to a suitable subgroup of Z*. This special structure of these solutions implies certain 
constraints on the possible range of values of the power spectral densities of the two sequences 
associated to the SDS. Specifically, we state and prove the fact that the power spectral densities 
remain constant over the orbits. 

Let (X,Y) be an SDS of Z v with parameters (v;r,s;X), with v odd and X = r + s — ^y-, 
corresponding to a circulant D-optimal matrix. In particular, we must have \X\ = r and 
\Y\ = s. Assume that 

X={JO„ Y={JO k (9) 

jeJ k&K 

for some subsets J, K of {1, 2, ... , m}. 

By abuse of notation, let X also denote the sequence x , xi, . . . , x v _i where 

f 1 if i^I 
Xi ~ \ -1 if ieX 
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and define similarly the sequence Y = y , yi, . . . , y v -\. 

Theorem 2 If k and k' belong to the same orbit O r C Z v and the sequence X is as defined 
above, then 

PSD x {k) = PSD x (k'). 



Proof 

Since k, k' belong to the same orbit O r we have k' = ks for some s G H . Since the set X is a 
union of orbits, we also have x as = x a for all a E Z v and s G H . Hence we obtain 



PSD x {k') = \DFT x {k 



/\|2 



v-1 

E 

a=0 



x a u 



tk' 



'v-1 

E 

v a=0 



X a UJ 




because x a G { — 1, +1} and lo = uj 1 . By expanding we have 

v-1 



PSD x {k') 



(a-P)k' 



v-1 

a =0,/3=0 
v-1 

^2 x a' > .' - 

a'=0,^'=0 



a=0,/3=0 
■u-1 



(Q-^)fcs 



u-1 

X as X/3s W 

a=0,/3=0 



(a-/3)ks 



(a'-P')k 



E 



a'=0 



I ft 'W 



v'A: 



\DFT x {k)\ 2 = PSD x {k). 



Here we have used the fact that when a ranges through Z V) so does a' = as, when s G H. 



□ 



Remark 1 Theorem^ says that the PSD is constant on the orbits of H . In fact it can easily 
be shown that PSD is constant on the orbits of the possibly larger group H* = H U (— 1)H. 
Theorem\^has some important computational consequences when one is using the power spectral 
density property as a preprocessor to discard candidate sequences and SDSs from a search for 
D-optimal matrices. In particular, it is clear that one needs to compute only one PSD value 
from each orbit Oj for 2 < j < m. Therefore, one can reduce the number of PSD values 
calculated per SDS from to at most m — 1 . 

Example 1 We illustrate theorem^ with a specific solution for v = 131, a prime. The auto- 
morphism group G of the cyclic group Z\^\ is cyclic of order 0(131) = 130. The subgroup H of 
G that was used was one of order 5. Its generator is 53, H = {1, 53, 58, 61, 89}. The H -orbits 
in Z v are {0} and the 26 cosets of H in G, thus m = 27. We enlarge this group by using the ele- 
ment — 1(= 130 in this case). Then we get a subgroup H* = {1, 42, 53, 58, 61, 70, 73, 78, 89, 130} 
of order 10. There are 13 cosets of H* in G, all of size 10 as well. Theorem^ implies that 

PSDx(k') = PSDx(k) and PSD Y (k') = PSD Y (k), W G H* ■ k. 

"'~ J — 13 PSD values for each sequence instead of = 65 



In this example, one can compute , 
values for each sequence. 
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4 The new circulant D-optimal matrices and SDSs 



We give D-optimal SDSs for the seven new parameter sets listed in the introduction. For some 
of them we give two or three non-equivalent SDSs. Non-equivalence of SDSs was established 
in all cases by an implementation of the method described in [4]. 

We define the notation for the orbits of the action of the subgroup that we use to construct 
the solutions below. Consider a fixed subgroup H of order h of Z*. Clearly h must divide 
\Z*\ = 4>{v). Denote by H ■ k the orbit of H in Z v through the point k, where • is multiplication 
mod v. We refer to the orbit H ■ = {0} as the trivial orbit. The orbit H ■ 1 is just the 
subgroup H itself. In general the size of an orbit may be any divisor of \H\ and if v is a prime 
then every nonzero orbit is just a coset of H in Z* and so will have size \H\. 
The notation 

X = |J H ■ j and Y = (J H ■ k (10) 

je.J keK 

will be used below to present all the solutions found, in fact each solution will be given only 
via the two indexing sets J and K. When X, Y are defined as above, then (X, Y) will be an 
SDS(t>; r, s; A) with r = J2 jeJ \H ■ j|, s — J2 keK \H ■ k\ and A = r + s - 

4.1 v = 63 

In this case the subgroup H is trivial and we write the blocks X and Y explicitly: 

X = {0, 1, 2, 3, 4, 6, 7, 11, 12, 13, 14, 20, 21, 22, 25, 26, 27, 30, 33, 35, 36, 38, 39, 42, 46, 48, 50, 53, 57}, 
Y = {0, 1, 3, 5, 7, 8, 10, 11, 13, 14, 16, 18, 19, 23, 30, 33, 34, 35, 39, 40, 48, 52, 54, 56}. 

This is an SDS with parameters (63; 29, 24; 22) which gives rise to a D-optimal matrix of order 
2 x 63 = 126. 

4.2 v = 93 

Consider the subgroup H — {1, 25, 67} of order 3, of Zg 3 . For the convenience of the reader we 
give below the 33 orbits of the action of H on Z g3 . 



H 


= 


= W 


H 


1 = 


= {1,25,67} 


H 


2 = 


= {2,41,50} 


H 


3 = 


= {3,15,75} 


H 


4 = 


= {4,7,82} 


H 


5 = 


= {5,32,56} 


H 


6 = 


= {6,30,57} 


H 


8 = 


= {8,14,71} 


H 


9 = 


= {9,39,45} 


H 


10 


= {10,19,64} 


H 


11 


= {11,86,89} 


H 


12 


= {12,21,60} 


H 


13 


= {13,34,46} 


H 


16 


= {16,28,49} 


H 


17 


= {17,23,53} 


H 


18 


= {18,78,90} 


H 


20 


= {20,35,38} 


H 


22 


= {22,79,85} 


H 


24 


= {24,27,42} 


H 


26 


= {26,68,92} 


H 


29 


= {29,74,83} 


H 


31 


= {31} 


H 


33 


= {33,72,81} 


H 


36 


= {36,63,87} 


H 


37 


= {37,61,88} 


H 


40 


= {40, 70, 76} 


H 


43 


= {43,52,91} 


H 


44 


= {44, 65, 77} 


H 


47 


= {47,59,80} 


H 


48 


= {48,54,84} 


H 


51 


= {51,66,69} 


H 


55 


= {55,58,73} 


H 


62 


= {62} 
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We give three non-equivalent SDSs with parameters (93; 45, 37; 36), via their index sets J, K to 
be used in ( flOl) . which give rise to D-optimal matrices of order 2 x 93 = 186. 

J = {2,5,8,9,10,12,13,24,33,36,37,40,43,47,55} 
K = {3,4,5,6,16,22,24,26,33,36,40,44,62} 

J = {3,6,8,10,11,12,13,17,18,22,29,33,37,43,55} 
K = {3,4,5,12,13,16,17,26,36,37,40,51,62} 

J = {3,5,8,10,11,12,13,16,17,20,24,29,33,48,55} 
K = {5,9,11,12,13,18,22,24,29,40,43,51,62} 

4.3 v = 103 

Consider the subgroup H = {1,46,56} of order 3, of Z* 03 . For the convenience of the reader 
we give below the 35 orbits of the action of H on Z W3 . 



H 


= 


= W 


H 


1 = 


= {1,46,56} 


H 


2 = 


= {2,9,92} 


H 


3 = 


= {3,35,65} 


H 


4 = 


= {4,18,81} 


H 


5 = 


= {5,24, 74} 


H 


6 = 


= {6,27,70} 


H 


7 = 


= {7,13,83} 


H 


8 = 


= {8,36,59} 


H 


10 


= {10,45,48} 


H 


11 


= {11,94,101} 


H 


12 


= {12,37,54} 


H 


14 


= {14,26,63} 


H 


15 


= {15,16,72} 


H 


17 


= {17,25,61} 


H 


19 


= {19,34,50} 


H 


20 


= {20,90,96} 


H 


21 


= {21,39,43} 


H 


22 


= {22,85,99} 


H 


23 


= {23,28,52} 


H 


29 


= {29,79,98} 


H 


30 


= {30,32,41} 


H 


31 


= {31,87,88} 


H 


33 


= {33,76,97} 


H 


38 


= {38,68,100} 


H 


40 


= {40,77,89} 


H 


42 


= {42,78,86} 


H 


44 


= {44,67,95} 


H 


47 


= {47,57, 102} 


H 


49 


= {49,66,91} 


H 


51 


= {51,75,80} 


H 


53 


= {53,69,84} 


H 


55 


= {55,58,93} 


H 


60 


= {60,64,82} 


H 


62 


= {62,71,73} 









We give three non-equivalent SDSs with parameters (103; 48, 42; 39), via their index sets J, K 
to be used in (fTUl) . which give rise to D-optimal matrices of order 2 x 103 = 206. 

J = {2,4,8,10,12,14,17,19,30,33,42,44,47,51,60,62} 
K = {1,2,3,14,20,21,30,33,38,40,42,44,53,60} 

J = {1,2,6,7,8,11,15,20,29,30,31,38,42,44,51,62} 
K = {1,4,11,12,19,21,23,29,31,38,40,47,51,53} 

J = {3,5,8,10,12,14,15,21,22,29,30,44,47,51,60,62} 
K = {1,2,5,14,15,19,23,29,31,33,40,42,51,55} 

We give one SDS with parameters (103; 46, 43; 38), via its index sets J,K to be used in (fTUl) . 
which gives rise to a D-optimal matrix of order 2 x 103 = 206. 

J = {0,2,7,10,11,12,15,17,19,29,31,33,38,40,42,47} 
K = {0,8,10,15,20,21,22,23,33,38,40,47,49,53,55} 



8 



4.4 v = 121 



Consider the subgroup H = {1, 3, 9, 27, 81} of order 5, of Z* 21 . We give one SDS with parameters 
(121; 55, 51; 46) which gives rise to a D-optimal matrix of order 2 x 121 = 242. 

J = {1,2,5,13,16,19,31,34,35,61,76} 
K = {0,7,8,10,13,16,22,25,26,40,76} 

4.5 v = 131 

Consider the subgroup H = {1,53,58,61,89} of order 5, of Z* 3l . We give two non-equivalent 
SDSs with parameters (131; 61, 55; 51), which give rise to D-optimal matrices of order 2 x 131 = 
262. 

J = {0,1,12,14,18,22,27,29,33,36,38,42,44} 
K = {2,4,8,11,12,17,22,33,36,38,42} 

J = {0,1,3,4,6,12,14,18,21,22,33,36,38} 
K = {2,3,11,12,17,18,19,27,29,38,42} 

4.6 v = 241 

Consider the subgroup H = {1, 15, 24, 54, 87, 91, 94, 98, 100, 119, 160, 183, 205, 225, 231} of order 
15, of Z£ 41 . We give one SDS with parameters (241; 120, 105; 105), which give rises to a D- 
optimal matrix of order 2 x 241 = 482. 

J = {3,4,5,6,7,10,13,38} 
K = {3,5,7,11,19,35,38} 

4.7 The Algorithm 

We now furnish a succinct description of the algorithm used in this paper. First we record the 
periodic autocorrelation function (PAF) values of the two candidate sequences A = [a , . . . , a v -±\ 
and B = [b , . . . , keeping in mind that the relationship 

PSD A (k) + PSD B (k) = 2v-2 

must hold for k — 1, . . . ,v — 1. This relationship is referred to as the PSD criterion and its 
importance lies in the fact that if for a certain value of k one of the PSD values is larger than 
2v — 2, then the corresponding sequence can be discarded, because of the non- negativity of the 
PSD values. The algorithm proceeds by generating randomly several millions of pairs (A, B) 
of ±1 sequences of length v that satisfy the PSD criterion. The A sequences are stored into a 
file and encoded by the concatenation of their PAF values on iJ*-orbits. The B sequences are 
stored into another file and encoded by the concatenation of their PAF values on iJ*-orbits, 
but where each PAF value is subtracted from the constant A. Then a distributed sorting algo- 
rithm is used to sort the two files in ascending order. Subsequently a linear time algorithm is 
used to detect any potential matches between the two sorted files. When a match is detected, 
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another algorithm is used to reconstruct the corresponding solution, based on some additional 
quantities that are computed together with the PAF values. 

We also mention here that in case v is a prime number, theorem [I] and its corollaries are not 
used in the algorithm, since they are not applicable. 

In terms of the algorithm performance, we have noticed that often it suffices to generate ran- 
domly anywhere between 10 to 100 million A and B sequences. The distributed sorting phase 
is quite fast, as it proceeds in breaking down the original files into several smaller files. The 
matching phase is also quite fast. The reconstruction phase is of negligible cost, computation- 
ally. 

In terms of the algorithm total running time for all D-optimal matrices found in this paper, we 
used approximately 200 CPU days to generate candidate A and B sequences and approximately 
an additional 100 CPU days to perform the distributed sorting and matching phases. 

5 Open cases for D-optimal SDSs with v < 200 

We summarize in the following table all odd integers v < 200 and the corresponding D-optimal 
SDS parameters, pointing out which cases remain open and which cases are known to exist. 
In constructing this table, we consulted information made available to us in [IB]. We also 
consulted and updated the table in [3J. 

Note that for v = 75 a non-circulant type D-optimal matrix of order 2 x v = 150 has been 
constructed in [TTJ. A circulant type D-optimal matrix of order 150 is not currently known. 
The asterisk in the "Existence" column indicates that such an SDS with the corresponding 
parameters is given in the previous section of this paper. 

Often there are more than one essentially different sets of solutions of the Diophantine equation 
a 2 + b 2 = 4v — 2, for a specific value of v. For each such set, we have elected to use the positive 
values of a, b, such that a < b, to construct the SDS parameters r, s, which implies that r > s, 
in all the SDS parameters that we present in our table. 
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